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A simple generalization of the Swift-Hohenberg equation is proposed as a model for the pattern- 
forming dynamics of a two-dimensional field with two unstable length scales. The equation is used 
to study the dynamics of surface waves in a fluid driven by a linear combination of two frequencies. 
The model exhibits steady-state solutions with two-, four-, six-, and twelve-fold symmetric patterns, 
similar to the periodic and quasiperiodic patterns observed in recent experiments. 
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Parametric excitations of surface waves have been ex- 
tensively studied since their first discovery by Faraday |l]] 
over a century and a half ago. In the basic experimental 
setup an open container of fluid is subjected to vertical 
sinusoidal oscillations, which periodically modulate the 
effective gravity. When the driving amplitude a exceeds 
a critical threshold a c a standing-wave instability occurs 
with temporal frequency u> one half that of the driving 
frequency. The characteristic spatial wavelength of the 
standing-wave pattern is selected through the dispersion 
relation uj(k) of the fluid. One typically observes patterns 
of stripes or squares in such experiments. It is only in 
recent years that a variety of additional patterns — some 
with quasiperiodic rather than periodic long range order 
- have been observed |§-§]< We shall focus here on a 
particular set of experiments, performed by Edwards and 
Fauve in which a fluid was driven by a linear combi- 
nation of two frequencies, forming periodic patterns with 
2-, 4-, and 6-fold symmetry, and quasiperiodic patterns 
with 12-fold symmetry. 

Previous theoretical work p|-pd|| has focused mainly on 
a description through amplitude equations with an angle- 
dependent interaction (3(0ij) between pairs of modes. 
Such an interaction, which is either postulated or derived 
from the underlying microscopic dynamics, can be chosen 
to stabilize iV-fold symmetric patterns for arbitrary N. 
Miillcr |l(]] has also used a set of two coupled partial dif- 
ferential equations, where the pattern of a primary field 
is stabilized by coupling to a secondary field which pro- 
vides an effective space-dependent forcing. Newell and 
Pomeau |nj have coupled multiple fields in a similar way. 
In both cases the coupling between the different fields is 
achieved through resonant triad interactions, similar to 
the interactions we shall introduce below. 

We propose a simple rotationally-invariant model- 
equation, governing the dynamics of a real field u(x, y, t), 
which describes the amplitude of the standing-wave pat- 
tern. Our approach is different in that it searches for the 
minimal requirements for reproducing the steady states, 
which are observed in the experiments of Edwards and 
Fauve We incorporate into our model only the two 
most essential aspects of the system: 

1. The dynamics is damped at frequencies away from 



the two forcing frequencies, and therefore the wave- 
lengths involved in the selected pattern lie in nar- 
row bands about two critical wavelengths. 

2. The driving used in the experiments is such that 
the up-down symmetry, taking u to —u, is bro- 
ken allowing interactions among triplets of standing 
plane-waves to exist. These triad interactions are 
the only stabilizing mechanism for non-trivial pat- 
terns in our rotationally-invariant model equation. 

We capture the essential dynamics with a single field and 
without a priori specifying any angle-dependent interac- 
tions among critical modes. This allows for a meaningful 
comparison of the stability of different A^-fold symmetric 
states. We find patterns of 2-, 4-, 6-, and 12-fold symme- 
try that are globally stable, but none with 8- or 10-fold 
symmetry, which is in agreement with the experimental 
observations of Edwards and Fauve || . 

The supercritical instability of a homogeneous state to 
a striped state is often modeled by the Swift-Hohenberg 
equation 12 
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which is variational, 
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driving the field u(x, y, t) towards a minimum of the Lya- 
punov functional (effective "free energy") — 
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The first term in the Lyapunov functional (g) favors the 
growth of the instability whereas the quartic term is re- 
sponsible for its saturation by providing a lower bound 
for T . The growth rate e of the instability is propor- 
tional to the reduced driving amplitude (a~ a c )/a c . The 
positive-definite gradient term is small only near the crit- 
ical wave number k c — 1, and thus inhibits the growth 
of any instabilities with wave numbers away from this 
value. 

If the parametric forcing is such that the u — > —u sym- 
metry is broken, then the Swift-Hohenberg "free energy" 
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is modified by the addition of a cubic term, — cra 3 /3. 
Such a term allows triad interactions of standing plane 
waves to lower the value of T and form hexagonal pat- 
terns. The analysis of the Swift-Hohenberg equation in 
the presence of this term is summarized, for example, in 
the review by Cross and Hohenberg p3| . With single- 
frequency forcing one cannot break the u — > —u sym- 
metry, but with certain combinations of two frequencies 
the up-down symmetry is broken and triad interactions 
become important. 

We model the two-frequency parametric excitation of 
a fluid by replacing the wavelength-selecting term in the 
Swift-Hohenberg equation (|l|) by a similar term which 
damps out all modes except those near one of two critical 
wavelengths: 
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The parameter c can be scaled out, but we include it here 
because it is used in the numerical simulations, shown 
later. Other model equations with similar wavelength- 
selection properties are possible. We choose this equa- 
tion because it is the simplest one that incorporates the 
physics we are interested in — it allows two unstable 
length scales and contains triad interactions among the 
different modes. Since (|j) can be applied to any pattern- 
forming system satisfying these requirements, it is not 
our intention to provide a detailed derivation of it from 
any specific underlying microscopic dynamics. 

Let us turn now to an analytic investigation of the 
model equation When both e and a are sufficiently 
small (or c sufficiently large) the wavelength selection by 
the gradient term is nearly perfect and the Lyapunov 
functional may be written in Fourier space as 



|k|=l,g 
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where the summations are restricted to wave vectors 
whose magnitude is either 1 or g, lying on two rings in 
Fourier space. The set of Fourier coefficients Uk, that 
give rise to the lowest value of T for a given choice of 
the parameters e and a, determines the most favorable 
steady state solution of the model equation ([!]). We are 
only interested in finding the global minimum of J-, thus 
establishing that our model indeed predicts the existence 
of the patterns observed in the two- frequency parametric 
forcing experiments. Of course, this approach may over- 
look meta-stable states or local minima of the free en- 
ergy. Note that with the omission of the gradient term, 
one may perform a rescaling of the field u — > cm. The 
rescaled free energy a~ A J- is then controlled by a single 
control parameter 
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FIG. 1. Phase diagram of the lowest-energy steady-state 
solutions of the model equation (Q) for q — 2cos(7r/12). The 
phase boundaries are lines of constant e* = e/a 2 . 



To study the formation of dodecagonal patterns we 
choose q = 2cos(7r/I2), which is the magnitude of the 
vector sum of two unit vectors separated by an angle 
of 30 degrees. We minimize the Lyapunov functional (|5|) 
with respect to the Fourier coefficients describing four 
different pattern candidates: (a) a striped pattern with 
space group P2mm, whose Fourier spectrum contains 
two opposite wave vectors of equal length; (b) a pat- 
tern of perfect hexagons with space group P6mm, whose 
Fourier spectrum contains a single 6-fold star of wave 
vectors; (c) a pattern of compressed hexagons with space 
group P2mm, whose Fourier spectrum contains four vec- 
tors on one ring and two vectors on the other ring; and (d) 
a dodecagonal pattern with space group P12mm, whose 
Fourier spectrum contains two 12-fold stars of wave vec- 
tors, one on each ring. 

We use standard methods |l4|] to calculate T for each 
of the cases. Because all the candidate patterns have 
symmorphic space groups [fl5[ which are also centro- 
symmetric we may always take all the Fourier coefficients 
on a given ring to be equal and their phases may all be 
chosen such that they are either or n. The minimiza- 
tion of the Lyapunov functional is therefore always with 
respect to no more than two real variables. We find the 
values of the Lyapunov functional for the different pat- 
terns to be 
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FIG. 2. Numerical solutions of the model equation (^) showing real-space patterns along with their Fourier spectra for 
different values of the control parameter e* = e/a 2 . The real-space images of u(x, y,t —* oo) show one quarter of the simulation 
cell with darker shades corresponding to larger values of the field. All figures are drawn to the same scale. In cases (a)-(d) 
q — 2cos(7r/12) = (2 + v / 3) 1//2 : (a) A 2-fold pattern of stripes for e* = 2; (b) A 2-fold pattern of compressed hexagons for 
e* — 0.1; (c) A 6-fold pattern of perfect hexagons for e* = 1.8; (d) A 12-fold pattern for e* = 0.015. In (e) q — 2, e* = 0.04, 
yielding a 2-fold superstructure of stripes. In (/) q = v2, e* = 0.04, giving rise to a 4-fold pattern of squares. 
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For e* > 1.91313 the striped pattern has the lowest free 
energy. For 1.91313 > e* > 0.08776 the 6-fold pat- 
tern of perfect hexagons and the 2-fold pattern of com- 
pressed hexagons (denoted by 4-2), which are degenerate, 
are most favorable. For e* < 0.08776 the dodecagonal 
pattern is the most stable. These analytical results are 



depicted in the phase diagram of Figure |l[ Note that 
the phase diagram depicts only the boundaries between 
global minima; in certain regions of the phase diagram 
additional states may be locally stable. 

The model equation (Q), supplemented with periodic 
boundary conditions, was solved numerically on a square 
domain using a pseudo-spectral method. The unit cell 
was typically chosen so that the simulation region held 
about 30 wavelengths. The simulation was performed 
on a 256x256 grid, with Adams-Bashforth second-order 
time-stepping. The value of c was taken to be be- 
tween 10 to 100. Figures H(a)-||(<i) show the real-space 
and Fourier-space results of the simulations with q = 
2cos(7r/12) for varying values of the control parame- 
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ter s*. The results are consistent with the Lyapunov- 
functional analysis and the phase diagram of Figure [|. 

Eight-fold and ten-fold symmetric patterns are not ob- 
served in our model for any choice of q. An analytic 
calculation of the Lyapunov functional (^|) for these pat- 
terns shows that it is greater than the free energy T§ 
© of the six-fold state, for any value of the control pa- 
rameter s*. This is in accord with the experiments of 
Edwards and Fauve ||, where such patterns are not ob- 
served. This does not rule out the possibility that octag- 
onal and decagonal patterns are locally stable but only 
that within the limits of our model they are not glob- 
ally stable. Two additional patterns that are observed in 
our model are a superposition of stripes of periodicities 
2ir and it (shown in Figure 0(e)) and a square pattern 
for q = \/2 (shown in Figure ||(f)). The latter has been 
reported by Edwards and Fauve. If one examines the 
Lyapunov functional (^J) in its full generality by allow- 
ing the value of q and all the amplitudes and phases to 
vary independently, other patterns might be discovered. 
We have only examined the symmetric patterns discussed 
here. 

The simplicity of our model shows that for continu- 
ous media very little is required to stabilize structures 
with quasiperiodic long range order: two length scales 
and triad interactions. The reason that 12-fold patterns 
are stable and 8- and 10- fold patterns are not is purely 
geometrical. In view of the Lyapunov functional , the 
crucial issue is the competition between the number of 
modes, which tends to increase the value of J 7 , and the 
number of triad interactions, which tends to decrease the 
value of T. The dodecagonal pattern of Figure |](g?) con- 
tains 24 non-zero Fourier modes and 32 distinct triangles. 
The octagonal and decagonal patterns do not contain a 
sufficient number of triangles to compete with the 6-fold 
pattern of Figure |](c). Our model confirms the conclu- 
sion of Edwards and Fauve that "12-fold patterns are 
more common than previously supposed." 

Our simplistic model is clearly not adequate for study- 
ing the structural stability quasicrystals in the solid state, 
yet it may offer a very simple system in which to study 
general questions regarding quasiperiodic order. These 
may include such questions as the formation and prop- 
agation of defects and phase boundaries as well as the 
dynamics of phason modes fl6|| . Moreover, we note that 
(01) may apply to situations other than Faraday waves. 
Any physical system that can be tuned such that two 
wavelengths undergo a simultaneous supercritical bifur- 
cation can be described by an equation similar to (^). 

An equation similar to (Q) could be used to study 
multiple-frequency forcing of Faraday waves with more 
than just two frequencies, as suggested by the title of this 
letter. We may speculate that with three or four forcing 
frequencies it might be possible to stabilize quasiperiodic 
patterns with even higher orders of symmetry, such as 18 
or 24. We leave the stability of higher-order symmetric 
patterns as an open theoretical and experimental ques- 



tion. 

We are grateful to Michael Cross, Jonathan Miller, and 
Peter Weichman for many helpful discussions. This re- 
search was supported by the California Institute of Tech- 
nology through its Division Research Fellowships in The- 
oretical Physics. 



[1] M. Faraday, Phil. Trans. Roy. Soc. Lond. 121 (1831) 299. 
[2] B. Christiansen, P. Alstr0m, and M. T. Levinsen, Phys. 

Rev. Lett. 68 (1992) 2157. 
[3] W. S. Edwards and S. Fauve, Phys. Rev. E 47 (1993) 

R788; J. Fluid Mech. 278 (1994) 123. 
[4] J. P. Gollub, Proc. Natl. Acad. Set. USA 92 (1995) 6705. 
[5] M. Torres, G. Pastor, I. Jimenez, and F. Montero De 

Espinoza, Chaos, Solitons, and Fractals, 5 (1995) 2089. 
[6] T. Besson, W. S. Edwards, and L. S. Tuckerman, Phys. 

Rev. E 54 (1996) 507. 
[7] B. A. Malomed, A. A. Nepomnyashchii, and M. I. Tri- 

belskii, Sov. Phys. JETP 69 (1989) 388. 
[8] W. Zhang and J. Vinals, Phys. Rev. E 53 (1996) R4283. 
[9] P. Chen and J. Vinals, "Amplitude equations and pattern 

selection in Faraday waves," preprint. 
[10] H. W. Miiller, Phys. Rev. E 49 (1994) 1273. 
[11] A. C. Newell and Y. Pomeau, J. Phys. A, 26 (1993) L429. 
[12] J. B. Swift and P. C. Hohenberg, Phys. Rev. A 15 (1977) 

319. 

[13] M. C. Cross and P. C. Hohenberg, Rev. Mod. Phys. 65 
(1993) 851. 

[14] See, for example, P. M. Chaikin and T. C. Lubensky, 
Principles of condensed matter physics, (Cambridge Uni- 
versity Press, Cambridge, 1995), ch. 4.7; or, for more de- 
tail, see L. Gronlund and N. D. Mermin, Phys. Rev. B 
38 (1988) 3699. 

[15] For a definition of "symmorphic space group" in the case 
of quasiperiodic structures and of the notation 12mm see, 
for example, D. S. Rokhsar, D. C. Wright, and N. D. Mer- 
min, Acta Cryst. A 44 (1988) 197. 

[16] H. W. Miiller, Phys. Rev. Lett. 71 (1993) 3287. 



4 



